Competition Between Pressure and Gravity Confinement in Lya Forest Observations 
Jane C. Charlton 1 ' 2 , Edwin E. Salpeter 3 , and Suzanne M. Linder 1 



Abstract 

A break in the distribution function of Lyct clouds (at a typical redshift of 2.5) has 
been reported by Petitjean et al. (1993). This feature is what would be expected from a 
transition between pressure confinement and gravity confinement (as predicted in Charlton, 
Salpeter, and Hogan (1993)). The column density at which the feature occurs has been 
used to determine the external confining pressure, ~ 10cm _3 K, which could be due to a 
hot, intergalactic medium. For models that provide a good fit to the data, the contribution 
of the gas in clouds to O is small. The specific shape of the distribution function at the 
transition (predicted by models to have a non-monotonic slope) can serve as a diagnostic 
of the distribution of dark matter around Lyct forest clouds, and the present data already 
eliminate certain models. 
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1. Introduction 

Clouds containing neutral hydrogen (with column densities ranging from 10 13 cm -2 
up to 10 22 cm -2 ) are observed through Lyct absorption in quasar spectra. In this Letter 
we restrict ourselves to the epoch with the most reliable data, redshift z ~ 2 to 3. The 
distribution function f(Njjobs) for observed HI column density of individual Lya absorp- 
tion lines is an important clue to the nature of the absorbing clouds, and the accuracy of 
determinations of / is improving steadily. With little loss of generality, an individual "Lya 
forest cloud" with some N Hobs can be modeled as a portion of a slab in hydrostatic equilib- 
rium with perpendicular HI column density Nh = Nnobs cos 6 (where 6 is the inclination 
angle) and an (unknown) combined HI and HII column density N tot (Charlton, Salpeter, 
& Hogan 1993 (CSH)). This modeling, in terms of a few adjustable parameters, is similar 
whether there are many different types of individual small clouds (minihalos) , each with a 
different characteristic value of Nh, or a single type of very extended disk structure (large 
proto-galaxies) with Nh decreasing enormously from the center outwards. 

Any realistic modeling of "forest clouds" must include the effects of both external 
pressure and self-gravity. CSH predicted a break in the distribution function, f (NHobs), 
where these two effects are equal. There must be a transition between pressure and gravity 
domination for some value Nniobs of column density. The most recent observational data 
(see Petitjean et al. 1993 (PWRCL)) shows a break of the type predicted, with Nniobs ~ 
10 16 cm -2 . The transition value predicted by the CSH models depends on the value of 
the external pressure, P ex t, and weakly on the ratio of dark matter to gas density at the 
transition; thus we shall be able to estimate P ex t- 

In this Letter we constrain model parameters by fitting to the neutral column density 
distribution over the forest regime. Since the observed distribution function f (NHobs) 
is of piece-wise power law form, it is a reasonable assumption that the basic distribution 
function g of the total (perpendicular) column density N to t is a continuous power law. The 
observed f(Nnobs) will be fit by adjusting three parameters: 1) the value of Nhi for the 
transition between pressure and gravity dominated regimes (set by the external pressure 
P ext and a dark matter parameter rj); 2) the slope of an assumed power law dependence 
of rj on N tot (The pressure induced by the self-gravity of the gas alone is proportional to 
A 2 ot , but the gravity due to dark matter present in the slab increases this pressure by 
a factor of rj.) ; 3) the slope of the power law g(N to t) distribution. We shall examine 
all of this three-dimensional parameter space, but give the specific example of extended 
proto-galactic disks in Section 3c. 

Alternative models for the observed "break" have been given by PWRCL in terms 
of four types of clouds and a transition from metal-poor to metal-rich systems, but these 
require more parameters than the CSH model. We also omit, in this paper, possible 
deviations from ionization and thermal equilibrium (Miralda-Escude and Rees 1993), since 
we wish to test the simplest models first. 

Section 2 uses the relationship between neutral and total column density in various 
regimes (pressure dominated, gravity dominated, Lyman limit, and damped) to schemati- 
cally illustrate the various features in the distribution function of Nh- Section 3 describes 
more detailed models that include corrections for orientation, extinction at the upper end 
of the forest regime, and consideration of the vertical slab structure. The best fit model 
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parameters (and plausible ranges) to the PWRCL column density distribution are pre- 
sented. Also, the number of systems expected in the Lyman limit regime and the total 
contribution of forest systems to the mass density of the Universe are computed. Section 4 
concludes with a summary of the constraints on the parameters. A value of P ex t is derived, 
and the implications of constraints on the dark matter distribution are discussed. 

2. Overview of the Different Regimes 

The distribution g(N tot ) is assumed to have a power law form. Written in logarithmic 
notation, g'(N tot ) = N tot g(N tot ), and we define the parameter e' by 

g'(N tot )d\nN tot = CN t ~J d\nN tot . (1) 

The distribution function f'(Nii) = Nh/(Nh) for neutral hydrogen, in logarithmic no- 
tation, can be derived if N H is given as a function of N tot . Since f'(N H )d\iiN H = 
g'(N tot )dlnN tot we have 

f'(N H )=g'(N tot )/k(N tot ), 

k(N tot )=dlnN H /d\nN tot . (2) 

To derive Nh as a function of N tot in the forest regime, we model Lya clouds as 
slabs in hydrostatic equilibrium and ionization/recombination balance. (For simplicity 
we neglect the presence of He.) The basic differences between behavior in the pressure 
and gravity dominated regimes can be illustrated by balancing the forces in a half slab 
approximation: 

^Gm H 2 vNl t + P ext = 2n tot kT (3) 

where n tot represents the total number density of H (at a "typical" vertical distance from 
slab center). On the LHS, the first term is the force per unit area due to gravity, and 
P ex t is due to the inward, external force. The outward force (RHS) is due to gas pressure 
within the slab, with the factor of two arising because slabs are highly ionized so that free 
electrons also contribute to this pressure. The parameter r\ is the enhancement factor for 
gravity over the gas column due to dark matter, taken to depend on N tot as 

v(N tot )= Vl (N tot /N tl )- 5 (4) 

where rji and N t i are values at the transition between the pressure and gravity dominated 
regimes (where the two terms on the LHS of equation (3) are equal). For the case with 
no dark matter we have r/i = 1 and 5 = 0. Using ionization/recombination balance for 
ionization rate ( (obtained from the proximity effect (Bajtlik, Duncan, & Ostriker 1988), a 
temperature of 20000K, and equating the gravity and pressure terms in equation (3) gives 
P ex t in terms of the observable Nhi(~ Nniobs/^-^) as 

- 1 ( 2 . 7xl C - 12s - 1 >< ife) ^ ™- 3 K. (5) 
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The relationship between N H and N tot in terms of their values at the transition, N H i and 
N tl , is 



where the subscript "hs" denotes the "half-slab" approximation used in this section. 

Of necessity there are four different regimes for the relationship between Nh and N tot 
in the range from the smallest to the largest values of column density, shown schematically 
in Figure la: As mentioned, there must be some value N t i of N tot (with a corresponding 
value Nhi of Nh) where the external and self-gravity pressures are equal. Far to the 
left of this transition (region (a) in Figure 1) the pressure is constant (= P ex t), Nh is 
proportional to N tot (eq. (6)) and k (eq. (2)) is approximated by the constant k a = 1. Far 
to the right of the transition (but still with Nh << 10 17 2 ) the ratio Njj/Ntot increases 
with increasing N to t as N^~ t & , so that in this gravity dominated regime (b) the derivative 
k is approximated by the constant kb = 3 — 5. The next transition point (2) occurs when 
Nh is near the value Nh2 = 10 17 ' 2 cm -2 where the optical depth for a (normally incident) 
photon just above the Lyman limit is 1. Just to the left of this transition, i.e. in the upper 
range of region (b), Nn/Ntot increases slightly more rapidly than N^ 5 because photon 
extinction decreases the ionizing flux and hence increases the neutral fraction. To the right 
of this transition, over a fairly narrow range of Ntot, from ~ Nt2 to Nt3 ~ 2Nt2, lies a 
region (c) which displays the "Stromgren effect", first suggested by Sunyaev (1969) and 
Felten and Bergeron (1969) and discussed in detail by Maloney (1993) and Corbelli and 
Salpeter (1993). The behavior of Nh in this region is complicated, but the essential point 
is that Nh increases by a very large factor, from Nh2 to Nh3 ~ N t 2 » Nh2- Finally, in 
the fourth region (d) where N tot » N tS: the gas is almost entirely neutral, Nh ~ N tot 
and kd ~ 1. 

If equation (1) held for all values of Ntot, the distribution function /'(Nh), given by 
equation (2), would also display the four regions shown in Figure lb. If we define j3' = (3 — 1 
(where f(N H ) oc N~ P ) by 



then well into regions (a) and (d), Nh oc N to t and k a = kd = 1 so that j3' a = j3' d = 
e' . Similarly, well into region (b), equation (8) holds with (3' b = e' /kb and kb = 3 — 
5. However, just before the transition from a to b the function k(N tot ) increases before 
settling into the constant kb, so that equation (2) predicts a short small dip, i.e. the 
slope din. f (Nh) / d\n.Nn does not decrease monotonically (see Figure 2a). The extinction 
correction (see section 3) near Nh2 has a similar effect of increasing k and there is again 
an increase in slope. In the Stromgren effect region (c) the extinction increases rapidly 
and N H /N to t increases rapidly at first and then approaches unity, so that k(N H ) increases 
at first, reaches a maximum and decreases towards unity. The actual behavior of f'(Nn) 
in region (c), indicated schematically by the dashed curve in Figure lb, depends on the 
shape of the photon spectrum and on just how large the ratio Nh^/Nh2 ~ N t 2/Nn2 » 1 
is (see Corbelli and Salpeter 1993), but f'(Njj) is not much smaller at the high Nh end 
of region (c) than at the beginning. We shall see that e' < 1 for any acceptable model, so 
that the total gas mass would diverge on the upper end of N to t if equation (1) held all the 




f'(N H ) oc N~ 



(8) 



4 



way up. There must thus be another transition point (4) beyond which g'(N tot ) changes 
and eventually decreases rapidly in region (e). 

The extent of the five regions (a) to (e) are not known a priori, nor are the transition 
values except for Nh2 ~ 10 17,2 cm -2 (defined by approximate unit optical depth at the 
Lyman edge). On the observational side there are three known transition points: Nh ~ 
10 13 cm -2 for the threshold of Lyo; absorption observations, Nh ~ Nh2 ~ 10 17 ' 2 cm -2 for 
the onset of the Lyman limit regime, and Nh ~ 10 20 cm -2 for the onset of the damped 
wing regime. Between the last two points, numbers of absorption lines can be measured, 
but not the actual values of Nh and we shall refer to this range from Nnobs = 10 17,68 cm -2 
to N Hobs = 10 20 - 5 cm~ 2 as "the long bin" of PWRCL. In principle, it is possible that 
Njji < 10 13 cm -2 so that all of the observable Lyct forest regime would lie in region (b), 
terminating at Njj2 ~ 10 17 ' 2 cm -2 . As Figure 1 shows, the predicted slope in this whole 
region (from 10 13 cm -2 to 10 17,2 cm -2 ) would have to be mostly constant with an increase 
at the end. The observational data of PWRCL suggests a decrease in slope, at best 
marginally compatible with a constant slope and incompatible with an increase in slope. 
We therefore accept their observational "break" near NHobs ~ 10 16 cm -2 as representing 
the theoretically required change in slope at N H i predicted by CSH. 

3. Fitting More Accurate Models 

We wish to fit all the observational data in the "Lyct forest regime" with log N Hobs 
up to 17.68, using eq. (2) but adopting several model refinements. 

a) Model Refinements 

(i) Instead of the "half slab approximation" of section 2, we solve for the vertical 
structure of a slab, i.e. the density p as a function of height z (from slab center). We 
assume throughout that dark matter gravity dominates over that due to the gas, and 
that the dark matter density is approximately constant over the height of the gas slab, 
so that the gravitational acceleration g(z) is proportional to z. With this assumption 
(plus assuming constant temperature) the shape of p(z) can be calculated in terms of the 
parameter N tot /N t i alone (see CSH, Section 5.1). 

(ii) In the approximate treatment so far, we have considered the ionization rate per 
H atom, a constant without allowance for absorption. By definition, observers chose 
for the upper bound to the "forest regime" (and lower limit to the "Lyman limit" regime, 
section 3d) a neutral column density Nnobs where the optical depth is moderate, r = 
3(13.6eV/-E) 3 for 10 17 ' 68 cm -2 ). In the small perturbation approximation, the percentage 
increase in Nn/Ntot equals the decrease in ( due to absorption, which depends only on 
Nh = N Ho bs cos 9 and the spectral shape (we use the usual quasar spectrum) . We have 
calculated this small correction (<^ 40%) for Nh up to 10 17 41 . 

{Hi) In reality, different slabs can have any inclination angle 6. The distribution of 
observed values of f (NHobs) can be derived from the distribution of perpendicular values 
(Milgrom 1988). Note that the orientation effect does not alter the power law behavior of 
f(N Ho bs) far from N H i (on either side), but is responsible for a gradual transition instead 
of a sharp break. 



5 



b) Fits to N H Distribution 

For any choice of the model parameters N H1 , 5 (eq. (4)), and e' (eq. (1)), the 
distribution /'(Nh) can be computed and normalized to provide the smallest value of x 2 
when compared to the PWRCL data. For various values of the dark matter parameter, 5, 
the best fit values of e' and Nhi are given in Table 1. The smallest x 2 value is obtained 
for 5 ~ —1.5, but the fit is only slightly worse for any 5 ^ 0.65. The slope of /'(Nh) as a 
function of Nh is illustrated in Figure 2a for various 5, and the f'(Nn) distributions are 
plotted along with the data in Figure 2b. Using g'(N tot ), which we have explicitly for the 
models in Table 1, the contribution to the cosmo logical density parameter is found to be 
Qf or ~ 10~ 2 1]^ 1 ^ 3 or slightly less for each of the four cases. 

c) Extended Proto-Galaxy Disks 

Hoffman et al. (1993) have proposed that the outer regions of proto-galaxy disks are 
responsible for the Lyman forest lines and that the total column density N tot (R) has a 
power law form. The exponent is related to e' in eq. (1) by 

N tot (R)~R-P, p=2/e'. (8) 

If we assume a spherical dark matter halo which gives a flat rotation curve, we can derive 
a unique relation between the parameters 5 and e': The gravitational force per unit area 
in the GDR (the first term in eq. (3)) is proportional to N tot W/R 2 where W = N tot /n tot 
is the slab width. Using equation (8) we find 

5 = 1-0.56', (9) 

Since e' in Table 1 is close to 0.71 (p=2.8) this gives 5 = 0.65. If the rotation velocity 
declines slowly at large distances 8 would be smaller. 

For assumed values for the constant rotation velocity V c and the radius -R13 where 
Nh has dropped to 10 13 cm -2 , one can calculate the transition dark matter parameter 771 . 
We find t]i = 4.7 (U c /70km/s) (_Ri3/250kpc) in terms of values for V c and R13 suggested 
by the "Vanishing Cheshire Cat model" (Salpeter 1993; Hoffman et al. 1993). For these 
values P ex t/k = 7cm _3 K, and the Lyct forest systems come from radii between 24kpc and 
i?i3, with a total gas mass of 8 x 1O 9 M per proto-galaxy. The dark matter mass (out to 
R13) is 3 x 1O 11 M per proto-galaxy, and the number density of these galaxies is about 
4.2Mpc -3 at z = 2.5 (the spheres would touch at z ~ 3.3). The contribution to O from 
the dark matter spheres is ~ 0.2. 

d) Model Predictions for Lyman Limit Systems 

In section 3b we obtained acceptable model fits to the observed distribution function 
f'(NHobs) up to log Nnobs = 17.68. In the "Lyman limit regime" we presently only have 
an observational value lib = 0.42 for the integral of / over a "long bin" stretching from 
log Nnobs = 17.68 to log Nuobs = 20.5 (PWRCL). Using equations (1) and (2) we calculate 
the integral of g' over the corresponding interval of N to t and equate it to the observed lib. 
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For each model in Table 1 we give the dimensionless ratio F[ b = Ii b e' / 'k(N t2 ) f (N H2 ) . If 
the power law in equation (1) holds for N tot up to N t3 , then is given by 

F lb = 1 - (N t3 /N t2 )- e ' . (10) 

The inequality Fi b < 1 (the requirement that the model predicts the observed number of 
Lyman limit systems for a finite value of N t3 ) is violated only for the model with S = 1.5 

IV. Discussion 

This letter has focused on fitting equilibrium slab models of Lya forest clouds to the 
neutral column density distribution of Petitjean et al. (1993) for clouds at z ~ 2.5. We 
find that the best fit is provided by models that have a transition between pressure and 
gravity confinement at neutral column densities of Nhi ~ 10 15 cm -2 . The best fit for this 
Nhi yields a x 2 °f 1-78 for e' = .70 and 5 = —1.5. However, fits are nearly as good for 
a larger range of parameters. Roughly, we have 14.8 < log Nhi < 15.5, .6 < e' < .8, 
and 5 < 1. In comparison, the best fit model in which the clouds are all gravity confined 
(N H1 < 10 13 cm~ 2 ) has e' = 1.3, 5 = .7 and % 2 = 3.5. Thus we maintain that the data 
requires a transition from pressure to gravity confinement. 

Since this transition is found to be at Nniobs ~ 2 x 10 15 cm~ 2 , equation (6) yields 

1/3 

Pext/k = 5.177]/ cm _3 K. This pressure could be due to a diffuse, universal, hot intergalac- 
tic medium. The measured value is smaller than limits placed by the COBE y parameter 
(distortion of the microwave background expected from Compton scattering by the IGM 
electrons) (Mather, et al. 1994). CSH favored a larger value of P ext based on extrapolating 
a single power law g'(N to t) to the damped wing systems, but region (e) in Figure 1 need 
not satisfy this power law. 

Models with 6 > 1 give a bad fit to the Lya forest systems and underpredict the 
number of Lyman limit systems (Fit, > !)• I n any case, the simplest models for the dark 
matter distribution suggest < 5 < 1 and the proto-galaxy disk model (with constant 
rotation velocity and a spherical halo) gives 5 ~ 0.6. However, a disk galaxy with a 
declining rotation velocity and/or a compressed spheroid for the inner dark matter halo 
would have an appreciably smaller value for 6. 

Although we have considered only simple models in this Letter, we hope that two 
features of our results are fairly universal: 1) The external pressure P ext is likely to be 
large enough to lead to a pressure/gravity transition at log Nhi ~ 15 and the neutral 
fraction of the hydrogen gas is likely to be large enough so that the gas in Lya forest 
systems contributes rather little to the cosmological mass density. 2) Because of the form 
of equation (2), the slope of the distribution function /'(Nh) should be a non-monotonic 
function (see Figure 2a). The present observational data in Figure 2b already seem to 
corroborate this feature. 

We acknowledge support from NASA grant NAGW-3571 at Penn State and NSF grant 
AST 91-19475 at Cornell, the careful reading of the referee, Phil Maloney, and the Aspen 
Center for Physics. 
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Table 1: Summary of Model Results 



5 e' logN H1 x 2 F lb 

-1.5 -70i;}o 15.0±;| 1.78 .21 

.7li;S)§ 15.1±;| 1.95 .37 

.65 .71 15.2+J 2.18 .57 

1.5 .69l;g? 15.4±;| 3.10 1.24 
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Figure Captions 



Figure 1: A schematic illustration of the regimes of behavior expected for Lya clouds. 
The relationship between the total column density N tot and the neutral column density 
Nh is given in the top half of the figure, and the lower portion sketches the resulting 
distribution function (the total neutral column density contained in systems with a given 
Nh)- The five regimes are: a) pressure dominated, b) gravity dominated, c) "Stromgren 
sphere" transition between optically thin and optically thick clouds, d) optically thick 
regime in which clouds are almost neutral, and e) high column density cut-off. 

Figure 2: The distribution function for the best fit models. The solid line represents 
the extended proto-galaxy disk model, a) The slope of the distribution function illustrates 
a feature in f'iNn) that is most pronounced for large negative values of 5. b) The distri- 
bution function is shown along with the data (la error bars) from PWRCL for the overall 
best fit model with 5 = —1.5, and for the proto-galaxy disk model. 
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